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^\ Abstract 

Chains of interacting non-Abelian anyons with local interactions invariant under the 
i— —i action of the Drinfeld double of the dihedral group D 3 are constructed. Formulated as a 

spin chain the Hamiltonians are generated from commuting transfer matrices of an integrable 
vertex model for periodic and braided as well as open boundaries. A different anyonic model 
with the same local Hamiltonian is obtained within the fusion path formulation. This model 

I is shown to be related to an integrable fusion interaction round the face model. Bulk and 

surface properties of the anyon chain are computed from the Bethe equations for the spin 
chain. The low energy effective theories and operator content of the models (in both the 
spin chain and fusion path formulation) are identified from analytical and numerical studies 
of the finite size spectra. For all boundary conditions considered the continuum theory is 
found to be a product of two conformal field theories. Depending on the coupling constants 

I the factors can be a Z<± parafermion or aMj 5 @j minimal model. 
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1 Introduction 

In recent years there has been a surge of attention directed towards the understanding of 
many-particle systems exhibiting topological order, i.e. phases which cannot be characterised 
by a local order parameter. Possible realisations of such topological quantum liquids in con- 
densed matter physics are the fractional quantum Hall (FQH) states [391 05] and certain two- 
dimensional frustrated quantum magnets (5j l35l 01] . The excitations in these systems display 
anyonic statistics and understanding their collective behaviour is essential for the classifica- 
tion of topological phase transitions. Particularly interesting are non-Abelian anyons where 
the interchange of two particles is described by non-trivial representations of the braid group 
complemented by fusion rules for the decomposition of product states. The fact that these 
non-Abelian anyons are protected by their topological charge has led to proposals for the use 
of such systems in universal quantum computation [361 US] ■ 

Some insight into the peculiar properties of many interacting anyons can be obtained in the 
context of simple model systems: such models can be obtained by associating anyonic degrees 
of freedom with each site of a lattice and defining interactions compatible with their braiding 
and fusion rules [T7]. The phase diagram of the resulting lattice models can be studied based 
on the numerical computation of finite size spectra. This approach is particularly powerful for 
anyonic chains, i.e. one-dimensional lattices, where the numerical data can be compared against 
predictions from conformal field theory (CFT). Another approach, also in one dimension, makes 
use of the fact that the lattice model may become integrable for particular choices of coupling 
constants. For the solution of such models various analytical methods, e.g. in the framework of 
the Quantum Inverse Scattering Method (QISM), have been established which allow to study 
the spectrum of their low-energy excitations, their thermodynamical properties including the 
long-distance asymptotics of correlation functions and even form factors |37} I38j. 

So far much of the work on such lattice models has been focussed on systems of the non- 
Abelian Ising or Fibonacci anyons related to the quasiparticles in certain FQH states and their 
generalisations appearing in su(2)k Chern-Simons theories [H [TTJ [281 113 ESS [5l] ■ These anyons 
have relatively simple fusion rules which allows for tractable computation of systems with 
nearest and next-nearest neighbour interactions. Integrable points within the one-dimensional 
versions of these models have been identified with restricted solid-on-solid (RSOS) or interaction 
round the face (FRF) models constructed from representations of Temperley-Lieb algebras |17| 
[3"3"j 132] . An alternative method to define an anyonic theory is via the Drinfeld doubles of a 
finite group algebra, D(G), and its representations [12]. The quasi-particles in these systems 
are irreducible representations (irreps) of D(G) labelled by their flux, i.e. an element of h G G, 
and their topological charge determined by the transformation properties under the residual 
global symmetry commuting with the flux h. 

Being a quasi-triangular Hopf algebra the quantum double allows for a direct construction of 
integrable quantum chains with nearest neighbour interactions described by a local Hamiltonian 
which is invariant under the corresponding symmetry I18j: within the QISM one obtains 
quantum spin chains on a Hilbert space being a tensor product of the finite-dimensional local 
spaces corresponding to a spin S, a qudit or a more general n-state quantum system. On the 
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other hand, it is already known that for any given model whose local Hamiltonian has the 
symmetry of a quasi-triangular Hopf algebra associated with an anyonic theory, it is possible 
to construct quantum chains using the fusion path formalism [19J. Here the basis vectors 
are composed of sequences of anyons and we shall refer to this as a fusion path chain. The 
local Hamiltonian is formally identical in the spin and the fusion path formalism. Therefore, 
one should expect the bulk properties of the spin and the fusion path model to be the same. 
The finite size spectrum of low energy excitations, however, is known to depend on boundary 
conditions [31 [8] and therefore should differ between the two realisations. 

In this paper we study this problem for a specific one-dimensional anyon chain with nearest- 
neighbour interactions. The underlying symmetry of the Hamiltonian is that of the Drinfeld 
double of a dihedral group, specifically D(D^). In the following section we define this algebra 
and recall its irreps and the corresponding fusion rules. Then, using the spin basis, integrable 
models are constructed subject to periodic, braided and open boundary conditions, all of which 
being based on the usual QISM transfer matrix. In the fusion path basis we use the approach 
of the previous work and find a fusion IRF transfer matrix whose series expansion contains the 
global one-dimensional Hamiltonian. In Section 3 we compute the bulk and surface properties 
of the model from the Bethe ansatz formulation of the spectral properties for the spin chain. 
The conformal field theory and operator content for the periodic spin chain version of the D(D^) 
model has been identified previously [20]. In Section 4 we expand this work providing more 
details on the analysis as well as extending the study of the finite-size spectrum to the spin 
chain with braided and open boundary conditions. In addition we present numerical results 
for the fusion path chain in support of the expectation that the low energy excitations of the 
D(D 3) -any on chain are described by the same CFT for all types of boundary conditions studied 
here. 

2 The Model and its Symmetries 

2.1 The D(D 3 ) algebra 

The model we consider in this article has the underlying symmetry of the Drinfeld (or quantum) 
double of a finite group algebra. The finite group we utilise is the dihedral group of order six, 
-D3, and is isomorphic to group of permutations on three elements, S3. This group is based 
upon the symmetries of an equilateral triangle and has the presentation, 

D3 = {<j,t\<j^ = = UTUT = e}, 

where e is the identity element of the group, a is a rotation and r is a flip. The Drinfeld double 
of this group is defined as the vector space, 

D(D 3 ) = C{gh*\g,heD 3 }, 

where * denotes an element from the dual space of CD3. This space forms a quasi-triangular 
Hopf algebra when equipped with the multiplication and coproduct, 

gihlM = S^l] (g l92 )h* 2 and A(gh*) = £ g^ti)* gk*. 

keD 3 

The remaining structure is uniquely determined by these relations [9j HI] . This algebra has an 
associated universal i?-matrix i.e. an algebraic solution to the Yang-Baxter equation. 
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Representations 

The representation theory of the Drinfeld doubles of finite group algebras are well known |14|l30j. 
The irreducible representations of D{D^) are classified by the conjugacy classes of D3. The 
irreps associated with the conjugacy class {e} are: 

tt^V) = 1 tt^Ct) = ±1 vr^Gf) = S e g 

,«-v> - ( : ) ^>w - ( ; 1 ) -(«•). 

where to = e~3~ . The irreps associated with the conjugacy class {a, a 2 } are: 

where 6 {0,1,2}. The irreps associated with the conjugacy class {r, err, a 2 r} are: 

/ \ / 1 \ / s T g \ 

vr^ ± )( C j) = 7 r( r ' ± )(r) = ±001 Tr^V) = Sf T 

V / V 1 / \ o o s° T ) 

The anyonic theory corresponding with D(Ds) associates an irrep with an anyon [12]. For 
convenience it is simpler to denote each irreps by a single letter, a, f). We equate 

a = (e, +), b = (e, -), c = (e, 1), = (a, 0), c = (a, 1), f = (a, 2), = (r, +), f) = (r, -) 

Properties of the anyons are inherited from their associated irreps, e.g. the dimension of an 
anyon equals the dimension of its corresponding irrep. 

Fusion Rules 

Required for an anyonic theory are the fusion rules of particles. These rules are defined by the 
tensor product decompositions of the associated irreps: 
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2.2 Local Spin Hamiltonians 

The D(Ds) model is constructed by taking a special case of the three state Fateev-Zamolodchikov 
model [16]. This limit yields the R- matrix, which can also be constructed from the 7r g ® 7r g 
representation of D{D^) [IB] . 

2 

R(z u z 2 ) = N( Zl ,z 2 ) [^ (i_i)(a " 6) ^(^i|a)^(^>)] e i+a+M® e y+a+6,i ; (1) 

a,6,i,j'=0 
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where e^j represents a 3 x 3 matrix (whose indices are considered modulo three) with an one 
in the ith row and jth column and zeros elsewhere, 



W(z\l) 



z-1 



and N(z±, z 2 ) = — -{wz\ — w 2 )(w — w 2 z 2 ). 



wz — ur 

The i?-matrix satisfies a Yang-Baxter equation in both spectral parameters, 

^12(^1,^2)^23(^12/1,^22/2)^12(^/1,2/2) = #23(2/1, yi)Ri2{xiyi, x 2 y 2 )R 2 z{xi, x 2 ), (2) 

and has the symmetry of D{D%), implying that the operator can be expressed in terms of 
projection operators. The projection operators from the 7T g (g>7r to the irreps in its decomposition 
are, 



V 



(a) _ 



9,h 



where tr is the trace. In terms of these projection operators the i?-matrix is written as, 

R(zi,z 2 ) = f a (zi,z 2 )pW + f c (z 1 ,z 2 ) p M +f X) (z 1 ,z 2 )p^ +f t (z 1 ,z 2 )p<- t) +MZ U Z2) P ®, (3) 
where 



fa(z 1 ,Z 2 )pW=R(z 1 ,Z 2 )pW, 



a e {a, c,9,e,f}. 



This R- matrix will allow the construction of integrable models subject to various boundary 
conditions. Due to the two spectral parameters present two local Hamiltonians are obtained in 
the usual manner, 

h {k) = i^-R(z u z 2 )-(3 k I®I, 

where k G {1,2} and G C is chosen such that the trace of the local Hamiltonians is zero. We 
find that, 



h (l) = 2^ p (a) _ V5 p (e) _ V3 p (5) _ 



2y/3 



3 

h m = 2 -fp^-fp 
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It follows that the local Hamiltonians commute with each other and with the action of the 
algebra: 



and 



{Tr s ®ir 3 )A(a),h^ 



0, 



for all a € D(D^). Therefore they have the underlying symmetry of D(D$) as the ii-matrix 
did. In this basis we have that and are related, 



where n is usual two-site permutation operator and * is, and herein reserved for, complex 
conjugation. Both local Hamiltonians are self-adjoint. 
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2.3 Global Hamiltonian 



We consider a variety of different models which differ only by boundary conditions. For every 
case the global Hamiltonian is comprised of two components in the following way, 



U e = cos(0) n w + sin(0) H 



(2) 



(5) 



The presence of the coupling parameter is due to the two distinct local Hamiltonians. For all 
boundary conditions considered the integrability of the global Hamiltonians is ensured by the 
existence of a commuting transfer matrix. Furthermore the two components of each global 
Hamiltonian will commute, 



o. 



(6) 



This commutativity allows us to investigate the individual spectra of T~L^ and 7i^ 2 \ Typically 
the spectra of and T~i^ will be identical or of a related form. 



2.3.1 Periodic Spin Chain 

We begin by considering the D(D^) model as a spin chain with periodic boundary conditions: 
its global Hamiltonians are defined by 

c-i 

n - n co + i(i+i)' 
i=i 

for k G {1,2}. Note that the periodic closure by the term in the global Hamiltonian breaks 
the D(Dz) invariance of the model. Both of these Hamiltonians appear in the series expansion 
of the commuting transfer matrix [38J, 



t(zx,z 2 ) 



tr [n £-Ro£(^i, z%) ■ ■ ■ IIoi-Roi^l, zjD] • 



By construction this transfer matrix is a polynomial of degree C in the variables z\ and z\ . It 
has been observed that this transfer matrix factorises and that its eigenvalues are always of the 
form 



A(zi, z 2 ) 



c c 
c \\[zi - z 1>e ) II O2 - z 2) e) 



(7) 



=1 



=1 



Therefore the eigenvalues can be conveniently described in terms of their zeroes Zk e = iuje Xk ' e , 
for k = 1, 2 and £ = 1, . . . , L. The two sets of parameters {xk,i}f=i independently satisfy the 
Bethe equations [21] 



(-1)' 



c+1 (l + (i/u)e*K 



1 — iuj e Xk <i 



n 

i=i 



to e Xk 'j 



j = !,...,£. 



(8) 



The energy eigenvalue of T~L^ corresponding to the set of Bethe roots {xk,i} is given by 

c 



EM = E({x k/ }) 



V 



1 



ioje Xk ' t 6 



3 + i\/3) C 



(9) 
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where we have used property sets of Bethe roots are invariant under complex conjugations, 
{ x k,t}f=i = {^fc^}f=i ( m °d 2wr). It is important note that while there are exactly C Bethe 
roots in each set {xk,e}f=i, they are allowed to be at ±00, but at most one at each. These roots 
need to be recorded and dealt with appropriately to ensure the energy is real. 

Let us remark that the Bethe equations Q and the corresponding energies ^ of the D{D 3 ) 
spin chain of even length C coincide with those of the three-state Potts spin chain with C/2 sites 
[2]. We shall use this equivalence below to identify some of the thermodynamical properties of 
the D(Ds) chain. 

The energy eigenvalues of the complete Hamiltonian are characterised by two solutions to 
the Bethe equations ([8]). As a consequence of Equations (5j6) and along with ^ they are given 
by 

E = cos(#) + sin(0) £ (2) = cos(0) E({x lyl }) + sin(0) E({x 2 , £ }), (10) 

provided these energies pair. Specifically, levels are said to pair if the two corresponding sets 
of Bethe roots form an eigenvalue of the transfer matrix, see Equation ([7]). As the two sets of 
Bethe roots need not correspond to a unique eigenvalue of the transfer matrix, e.g. there may 
be two eigenvalues that differ by a constant factor or an eigenvalue might be degenerate, we 
refer to the total number of eigenvalues, including degeneracies, as the pairing multiplicity. 

The total momentum of the corresponding state can also be given in terms of the two sets of 
Bethe roots: at z\ = 1 = Zi the transfer matrix becomes a shift operator by one site, therefore 
the momentum operator is P = — i lnt(l, 1). Unlike for the energy ^ the partial momentum 
contribution from one of the participating Bethe configurations cannot be defined uniquely [21] . 
By construction the eigenvalues of P are real (2tt/£ times an integer for periodic boundary 
conditions considered here) . Therefore we may define the partial and complete momenta to be 



P( fc ) = P({z M }) = Re 
P = pW _ p(2) + const _ _ 



iwe Xk ' e ) 



(11) 



Again, roots x^/ = ±00 have to be taken into account to ensure finite (partial) momentum. 
The total momentum is given by the difference of partial momenta reflecting the fact T~L^ is 
the spatial inversion ofHW. The remaining constant represents a macroscopic effect, details 
of which have been discussed in earlier works |21|. 



2.3.2 Braided Chain 

One closed chain proposed as an alternative to the periodic chain is the braided chain [22" l l3T | IM]. 
In this model, translational invariance is replaced by invariance under a global braiding operator. 
As a consequence the underlying symmetry of the model will not be broken, i.e. it has the full 
global D(Dz) symmetry. The global Hamiltonians are defined by, 



c-i 



fee {1,2}, 



3=1 



where 



B = fol2^23---^(£-l)£) 



and b = lim 

z— >oo 



;R{Z,Z) 
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There are different possible definitions for the braiding operator 6j relating to other limits of 
R(zi, Z2) which we do not consider within this article. The integrability of the model is ensured 
by the existence of a transfer matrix t(zi,Z2), which can be found in |18| . Eigenstates of this 
model are again characterised by the Bethe Equations Q. Again there must be C Bethe roots, 
this time with one Bethe root allowed at +00 but none allowed at —00. In this case we find 
t(l, 1) = B. Calculating the order of B for small system sizes we find that 

B n = j*l when ( n = 3£, C even, 
" n = 2£, £odd. 



This allows us to define a generator of translation, P = —ilnB. The eigenvalues of P are 

3C 



restricted to integer multiples of either j or |^ depending on the parity of C. 



2.3.3 Open Boundary Conditions 

We also consider spin chains with open boundary conditions. The reflection matrices, i.e. 
_RT-matrices, required for the construction of an open chain whose boundary fields preserve 
integrability [51] are the same as have been determined for the D(D^) one-parameter i?-matrix 
|10| . The global Hamiltonians are, 

^=X k B?~ +X tB^ + Y.^- 

The boundary operators B^~ (B^ + ) act on the first (last) site of the chain, respectively. 
There are three possible (independent) options for each of the different boundary operators, 



w 2 b w 2 b 2 
B^-,(BW + y,(B^-y,B^ + g { [ wb 2 b 

wb b 2 



b = 1, w, w 2 



The real boundary amplitudes where Xh > ^ e {1>2}, have to satisfy X1X2 = X1X2 = 0. Like 
the periodic and braided models integrability is derived from the existence of a transfer matrix 
(see |18t [2T] for the open D(D^) transfer matrix) and the eigenstates of the Hamiltonian are 
classified by sets of Bethe roots. The Bethe equations for the Hamiltonian T~L^ are independent 
of the choice of the boundary operators, 



yr ( e Xk > 1 - u 2 e Xk > A _ c+1 ( 1 + we 21 ^' Ul- uj 2 e 2xk ^ Ul + iw 2 e*w \ 2C 

11 I p^ fc ,i - ue xkj j ~ \ > V 1 + uj 2 e 2xk <i Ml- oje 2xk >i 111- ioje Xk '- 



1 = 1 



x ®(xk,j,X k )®( x k,j,Xk) 
with Bethe roots always appearing in pairs of ±x and where 



®( x k,j,Xk) = { ( X±v^(l-^ 2 ^)+^ 2 (l+xJV3K fe .i 
X ^V3(l-^e 2Xk ^)-iu J (l+x : tV3)e Xk ^ 



ii xt = 0, 



As the Bethe equations depend on the parameters so do the number of Bethe roots. If Xk = 
then we obtain the open chains with free-ends, in which case we have dk = 2C Bethe roots with 
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at most one pair of roots at ±00. We refer to this Hamiltonian as the bulk Hamiltonian and 
describe the Bethe roots of the system as belonging to the bulk. If only one of either xt 01 
is non-zero then there are df. = 2C + 2 Bethe roots, while if both \t an d Xu are non-zero 
then there are dk = 2£ + 4 Bethe roots. The extra pairs of Bethe roots appearing, as compared 
to the free-ends case, are referred to as boundary Bethe root^] The boundary Bethe roots 
are always finite but approach ±00 in pairs as Xk approaches zero. The energy of the 
Hamiltonian corresponding to a set of Bethe roots is given by 

£ (fc) (W» = ^{e 

1 1=1 

where 

Only when the boundary interactions are not present, i.e. with free-ends, will the model have 
D(Ds) invar iance. 



iwe 



2.4 Fusion Path Analogues 

As the local Hamiltonians have the symmetry of D(D%) it is possible to create fusion path 
analogues [19] . Depending on the boundary conditions imposed the global Hamiltonians may or 
may not be equivalent to their spin formalism counterparts. The construction of the analogous 
fusion path chains uses the same method of Pasquier's representation theory reliant face- vertex 
correspondence |49j . This allows the fusion path analogues to be considered as the Hamiltonian 
limits of RSOS models and proves their integrability. The connection between fusion IRF 
models and many other physical systems has already been established [25j . 

We first define the fusion path basis. Basis vectors of the fusion path space are of the form, 

\a ai...a c ) , 

where a, G {a, b, c, 0, e, f, g, b,} and neighbouring labels satisfy the condition, 

OiOi+i e {ab \VbCV a ® V g } (12) 
= {ag, bf), eg, cf), Dg, t)h, eg, ef), fg, fh, 
ga, gc, qD, ge, gf, fib, fjc, f)Z>, f)e, fjf}. 

Thus aj+i must appear in the fusion of aj and g. Diagrammatically a basis vector corresponds 
to the figure below, where the joining of two lines indicates fusion which occurs from top-left 
to bottom-right, 



a 



_04_ 



a 2 



«3 



ac-i 



l It is important to note that for \k sma ll the presences of the boundary Bethe roots will not have a significant 
effect on the configuration of the bulk Bethe roots. 
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To construct operators on this space we utilise F- moves (generalised 6-j symbols), which allow 
the temporary re-ording of fusion, 

be b 



5>< 



abc\d 
e Id 1 a 




These F-moves are required to construct the projection operators acting on a boncj^J 

[Fa i+1 V [F ai+1 ) b ^.ai-iOiOi+i..) (..ai-iaiai+i. 



v {b) 



a i-l, a i> a ,i a i + l 



where the coefficients F are the unitary F-moves associated with D{D^). These can be calcu- 
lated explicitly from the representation theory of D(D^), however, as they are not essential we 
omit them. Below we shall consider local interactions defined analogously to Equation Q: 



2\/3 

2y/3 



(a) 



-P 



(a) 



3 F 
V3 



(0 



\/3 \/3 _m , 2\/3 

~3~ P 



p(e) + 



(f) 



3" 3" 3" 3" 3 

These Hamiltonians act on three consecutive labels but only can change the middle label. 

As a consequence of the equivalence of the local interactions between the spin and fusion path 
formalisms the global models in the two formalisms may differ only by boundary conditions. 
The open model with free ends and braided model both have D(D^) invariance which means 
that the fusion path and spin versions of these chains are equivalent. For periodic boundary 
conditions neither the spin nor fusion path formalism has the complete D(D^) invariance and 
thus the two models, while sharing bulk properties, are distinct [19]. For a periodic model in 
the fusion path basis we need to consider the space spanned by the basis vectors satisfying 
a o = a>C- As a consequence of the fusion rules this periodic closure is possible for C even only 
and leads to the decomposition of the Hilbert space, 

Hilbert Space = C {\a±. .. ac) \a± = g or a\ = h} © C{|oj... etc) 1^2 = S or °2 = > (13) 

where each of these subspaces has dimension 3^ + 1. The global Hamiltonians are, 

c-i 



= h 



(£-!)£! 



h 



(fe) 
£12 



(3-1)3(3+1)' 



k G {1,2}. 



The integrability of this model arises from the existence of an analogous ii-matrix connected 
to an RSOS model whose heights correspond to the labels of the irreps of D(D%). As the R- 
matrix of Equation ([3]) is expressible in terms of the D(D^) projection operators it follows that 
there exists an equivalent operator in the fusion path basis 

R(zi,z 2 ) 



fa(zi,z 2 )p {a) + fc(zi,z 2 )pU + Mz 1} z 2 )p^ + / t ( Zlj Z2 )pM + / f ( Zl) Z2 )p(f\u) 



E 



B 



ai,a,2,a 2 ,a3 



a 2 0,3 



zi,z 2 ) \aia' 2 a 3 ) (aia 2 a 3 | 



2 It is important for the reader to note that in this fusion path formalism that the labels in the basis vectors 
do not correspond to individual sites but rather bonds. The individual sites are still the g-anyons but now can 
not be solely acted on as this would break local D(Ds) invariance. 
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which satisfies a face Yang-Baxter equation, 

#123(^1, ^2)^234(^12/1, ^22/2)^123(2/1, 2/2) = -R234G/1, 2/2)^123(^12/1,^22/2)^234(^1,^2), 
as R(z±, Z2) satisfies Equation fl2J). The weights appearing in R(z\, 22) are used to construct the 



commuting transfer matrix [29], 

(a' x ...a' c \ t(z 1 ,z 2 ) \a\...ac) 







21,22 



Once again the global Hamiltonians is generated by this transfer matrix, with and 
appearing in its series expansion, implying integrability. The existence of the transfer matrix 
also guarantees commutativity of the two components 

n^\n^] = 0. 

In previous anyonic fusion path models |17| [33] integrability was also shown by constructing 
transfer matrices associated with RSOS models, however in these instances the fusion path 
i?-matrices corresponded to representations of the Temperley-Lieb algebra. 

Every RSOS model can be naturally associated with a graph [35]. For the D(D^) model 
considered here we obtain the graph given in Figure [TJ This graph is equivalent to McKay's 
representation graph for the representation 7r g of D(D^) [53]. The graph indicates that the 
D(Ds) fusion path model does not correspond to any of the known (dilute) RSOS models 
associated with Dynkin diagrams j35J H?l [5U1 G2] • We notes that the graph also doesn't appear 
amongst the more general graphs associated with other RSOS models |13j . 



Figure 1: A graphical representation of allowed neighbouring labels in the fusion path chain. 
The vertices/nodes of the graphs are the labels of anyons which are connected via an edge if 



and only if the two anyon labels can appear next to each as given by Equation (12) 




3 The Bethe Equations and Exact Results for Spin Chains 



As a consequence of Eq. ([10]) the ground state energy of the model is always obtained by the 
following combinations (NB: provided that these states are allowed to pair) 



En 



(l) 



, + sin(0).E z (2) 

(i) 
h 

cos(6)e£ ] +sm(9)E { ^, 
{ cos(e)E { l 1) +sm(9)Ef\ 



cos(e)E i 

cos(9)E { h 1) +sm(e)E^\ 



< 

£ < 

2 — 

vr < 



< - 
^ 2 ' 

<7T, 



< 



3tt 



(15) 



3- 



< 



2 ' 
< 2vr, 
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where E^ is the lowest energy of the T-L^ and E^' is the highest. An immediate implication 
of this form of the ground state energy is that there are level-crossings for 9 being integer 
multiples of 5 leading to first order quantum phase transitions. Here we will use a different 
consequence of (15): the complete spectrum of the model can be obtained from an analysis at 

An 



these particular points in combination with the implementation of the pairing rules 
additional simplification arises from the fact that the spectra of and are identical for 
most boundary conditions considered in this paper: this allows us to restrict the analysis of 
the low-energy spectrum to those of the Hamiltonians %e=Q = T~L^ and He=n = —H^ whose 
ground state energies are E- 



(i) 



E\ 2) and -E 



(i) 
h 



(2) 

E) , respectively. Only in the case of 



open boundary conditions with D(L>3)-symmetry breaking boundary fields are the spectra of 
±U^ and ±ft( 2 ) are independent. 



3.1 Energy Density in the Thermodynamical Limit 

In the thermodynamic limit C — > oo bulk properties of the system are independent of the 
boundary conditions imposed. Therefore we can compute the energy density from the Bethe 
equations Q for the periodic spin chain. To this end the solutions to the Bethe equations 
need to be classified and the root configurations corresponding to the ground state and low 
energy excitations have to the identified. As mentioned above, the Bethe equations for the 
periodic D(D^) spin chain arise also in the context of the 3-state Potts model. For the latter 
the classification of solutions has been obtained by Albertini et al. [2j H] ; see also Ref . |21j . In 
particular, numerical diagonalisation of the transfer matrix shows that the lowest energy states 
of 7-Le=iT consist of three different Bethe root (zi = zij) types: 

1. Positive Bethe roots (+- string), z\ = iuje x ~i , where xf G M, 

2. Negative Bethe roots ( — string), zi = iu)e x i +m , where x7 G K, 



3. 2-strings, where the Bethe roots come in pairs, z\ 



X s -— ■,, 

iue t 3 with 



xf G 



as well as a limited number of Bethe roots at ±oo. Letting N + , N- and A2 be the number of 
+- strings, — strings and 2-strings respectively and setting n± oc G {0, 1} to be the number of 
Bethe roots at ±00 then we have the constraint 



N+ + N_ + 2N S = C- n , 



In the 3-states Potts model these root types were also identified, along with a few other which 
we don't consider, with the additional constraint n +oc = n_oo [2]. 

Allowing combinations of these roots we then find that the for the periodic Hamiltonian we 
can take the logarithm of the Bethe equations and define the following set of counting functions: 



Z+(x) 
Z-(x) 
ZJx) 



7 1 N+ 

1=1 

1 1 N+ 



(X — X 



12' 

11. 



1=1 

(■<' - ■'•/' ; + 7 2v 



(X — X 



I ' 



1 Ns 2 



1=1 
N. 



1=1 
1, 



1 > 6 . 
1 N+ 



(X — X 



1=1 



1 '3 

N. 



1=1 



c 



x — x 



1 ' 



1 N ~ 1 1 Ns 



k=i 



C 



1=1 



[x Xj : 



1=1 



'3' 
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where 



l-i f tanh(f; 

(x; t) = tan ; — - 

^ ' vr V tan ( i7r ) 



We find that the ground state for T~Le=-w consists entirely of 2-strings and is only realised when 
C is even. 

Similarly we find that the lowest energy states of %e=o consist of the three same Bethe root 
types and hence we have the same counting functions. However, we find that the ground state 
this time is only realised when £ is a multiple of four and consists of 3C /4 negative Bethe roots 
and £/4 positive Bethe roots. 

Based on these observations the root density formalism [57] can be applied to compute the 
corresponding energy densities: The density of 2-strings in the thermodynamic ground state of 
He=tT and their dressed energies 7-Le=-w are determined by linear integral equations 



p(x) 



e(x 



1 



1 



1 



\4cosh(rc) 
1 



2^/3 2 cosh(x; 



+ 



V3 
2vr 



dy 



2 cosh(x — y) + 1 



p(y) , 



1 



+ 



dy- 



1 



(16) 



<y) 



4 cosh(x) - 2-v/3 2 cosh(x) 2vr 7_ oc 2 cosh(x - y) + 1 

These equations (and the corresponding ones for Hq = q) can be solved by Fourier transformation 
giving the ground state energy densities [21 EE] 



Eft 



1 2^/3 



and 



C 



1 

2^ 



2^/3 3 
~9~ + 4 



(17) 



3.2 Fermi- velocity 

The low-energy excitations over these ground states have a linear dispersion and their Fermi 
velocities have been computed from the root density formalism in the context of the three-state 
Potts model The identification of the energy eigenvalues of this model with the partial 
spectrum of the D(D^) Hamiltonian ^ based on the equivalence of the corresponding Bethe 
equations however can not be extended to the spectrum of the (partial) momentum operator. 
Instead we use the observation that the root configurations solving the Bethe equations are 



invariant under complex conjugation {x^} -R- (mod 2iir) to rewrite (11) as 

2i ^-^ \\ — {l/iuj)e Xe 



(18) 



which is half of the momentum in the three state Potts chain. This implies that the contribution 
of a 2-string to the partial momentum is 



/X 
dyp(y) 



and therefore the Fermi velocity of low lying excitations of Hq =7T is found to be 

1 e'(x) 



(19) 



(20) 



7T p{x) 

in agreement with the finite size analysis of the spectrum performed in [20| . 

Similarly we can compute the Fermi velocity of gapless excitations for the Hamiltonian Hq. 
Again the result is twice than what has been found for the three-state Potts chain [2], i.e. 

v F = l 
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3.3 Boundary Fields 



We can also determine the exact expressions for the surface energy, i.e. the C° contributions 
to the energy, for the open model with interacting boundary fields. Firstly we find the ground 
energy for the Hamiltonians %e=-K an d T-ie=o, and then extend this result to generic 9. We can 
use the Bethe equations given in the previous section and apply the same method that was used 
to calculate bulk energy density, however, due to the free parameters we need to utilise certain 
assumptions concerning the behaviour of the boundary Bethe roots. Due to the symmetry of 
the Bethe equations it is enough to look at the correction to the surface energy with only one 
boundary field present compared to the free-ends case. 

As stated earlier for Tig^ with \i = and xt 7^ there will be two additional Bethe roots 
present. If xt is n °t too large then we find that ground state consists of C 2-strings (for even 
C) bulk Bethe roots and two boundary Bethe roots as either i-strings depending on the sign 
of xt- The smaller xt is, the larger the magnitud^jof the boundary roots are. Using these 
observations we can obtain analytical expressions for the surface energy for — 3^/3 ^ ^fe ^ ^73' 
We find that the correction to surface energy for the Hamiltonian ~H.q =v with one interacting 
boundary field, compared to the free ends case, is 



18x 3 \/3 



2xV3 

1 + 2x\/3 tt(1 + 2 X \/3) y/l + 2x\/3 - 9 X 2 



arccosh(-| - ^=), X < 0, 
0, " ' x = o, 

arccosh(i + ^75), X > 0. 



Similarly, we find the surface energy correction for T~L 



6»=0> 



9e=o(x) 



1 + 2 X V3 



+ 



9 X 3 v^arccoshc- § - ) 



2(l+2 X v / 3)Vv / 3-3x ^(^x^) V^X^^x 2 ' 

o, x = o, 



+ 



9x 3 V3arccosh(i+ 5 ^_) 



2y / v / 3+9x 7r(l+2xv / 3)v / l+2xv / 3-9 X 2 ' 

From these we are able to determine the ground state energies of Hg =7r and T-Lq=q 



X<0, 

x = o, 
x>o. 



E e= n(xt>Xi 
E e=0 (xt,Xi 



1 2^3 
7T + ~~9~~ 



C + 



3 2^/3 
2 3~~ 



+ g e=n (xt) + ge=n(Xi ) + o(C G 



1 2y/3 3 
2^ 9~~ + 4 



C + 



3 2^3 



(21) 



+ ge=o(xt) + 9e=o{Xi) + o(C°). 



In Figure [2] we plot the predicted functions g e=w (x) and ge=o{x) compared to the equivalent 
numerical values obtained by solving the the Bethe equations (with X\ = X2 = 0)- 
It is clearly seen that the numerical values match the analytically predicted results even outside 
of the region — < x < -4=. We also note that increasing the strength of the boundary 
field such that it is the dominant term gives the energies associated with the eigenvalues of the 
boundary operators, 



lim 1 

x^+00 ' 



Ax) ~ 2 X 



lim 

X->—o< 



ge=o(x) 



and lim g e=1 r{x) ~ ~X ~ lim < 

X->-oo X->+°° 



=0 



(x). 



On a technical note, it was observed the magnitude of the boundary Bethe roots will also increase as C, 
increases for fixed \. We assume that as £ goes to 00 the boundary Bethe roots tend to ±00 for fixed boundary 
amplitudes \k ■ This is important when considering the thermodynamical limit. 



14 



Figure 2: Numerical values of the functions g e=n (x) an d ge=o(x) obtained via solving the Bethe 
Equations compared to their predicted analytical expressions for a system of 100 sites. In both 
plots we have data points for x < — ^75 ~ —0.193 and x > -75 ~ 0.577. 




These limits show the correct asymptotic behaviour of the ground state energy, Equations (21 ), 
with respect to the x- This leads us to conjecture that the energies are analytically correct for 
all values of Xk e ' 



Using Equation (15) along with the relations, 



= -E e=7r (xt, X k ) and E\ k) = E e=Q { X + k , X~ k ), 

we are able to determine the ground state energy for of Hq for generic 9. We should again recall 
that we have the constraint xfxt = X1X2 = an d that the interacting boundary terms break 
the D(D%) invariance of the model. 



4 Excitations and Conformal Field Theories 

The presence of a coupling parameter makes the identification of a conformal field theory a 
difficult takes. Every energy level of the Hamiltonian is dependent on the coupling parameter 
9, moreover the energies which are the lowest will change with 9. This implies that model will 
not be described by a single conformal field theory but rather but multiple ones. To simplify 
the issue we first study the model at level crossings which turn out to be described by minimal 
models. This allows us to use powerful machinery and well-known results to accurately describe 
the model at these points. 

To then recover the complete model once all the level crossing have been identified is a 
relatively straight-forward, albeit non-trivial, task. The main difficulties comes with identifying 
the previously mentioned pairing rules and connecting them with some conserved quantity of 
the model. Once all this has been achieved the lowest lying energies will be known for all values 
of 9. 
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4.1 Periodic Spin Chain 

Here we present an extended version of the results found in the letter [20] . At the level crossings 
conformal invariance predicts the scaling behaviour of the ground state energy should be 

E = eooC - (cv F ) x + o(£ _1 ), 
bL 

where c is the central charge. The operator content of a given realisation of the CFT is 
constrained by modular invariance of the partition function and boundary conditions [7]. 
Further constraints are imposed by locality of the physical fields which considered for the 
complete model, i.e. not at a level crossings. The difference in energy and partial momentum^ 
must be of the form, 

E(C)-E (C)= 2 ^(X + n + n) , P(£) - P (£) = ^ (s + n - n) (22) 

the scaling dimensions X = h + h and conformal spins s = h — h of the primary fields in the 
theory can be determined numerically (re, re are non-negative integers). To obtain the data 
numerically we can use solutions to the Bethe equations along with Equations ^ and (11). 



4.1.1 Spectrum of %g=7r 

It is known that the ground state energy of H n is given by [2] 



En 



1 2^/3 
7T + ~~9~~ 



£ "I X ^ + ° (rl) (23) 



Using the Fermi- velocity we have that the central charge of the effective field theory for the low 
energy degrees of freedom in Hg =n is c = 4/5: this sector of the model is in the universality 
class of the minimal model ^(5^), the conformal weights h, h of the primary fields can take 
the rational values from the Kac table 



h, h G 



(6p — 5q) 2 — 1 
120 



1 < q < p < 5 



1 1 1 2 21 2 7 13 

'40' 15' 8' 5' 40' 3' 5' ~8~' 



The excitations occurring in this sector for even chain lengths are given in Table [T] 
For the states given we have solved the Bethe equations up to a minimum of 40 sites, although 
in general over 100 sites were considered when possible. As periodicity breaks the D(D^) 
invariance of the model we can no longer use the decomposition 

n ®£ = + ^ e _( 3 £-2 _ ^ e 3 C-2 ^ e ^ e ^ e ^ ^ 

for C even. However we find by combining some of these irreps we have residual symmetry 
allowing us to partition the excitations. Specifically the irreps tt c and irt can not be separated. 



At the level crossing there is massive degeneracy and the complete momenta are not unique but partial ones 

are. 
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Table 1: Scaling dimensions X v extrapolated from the finite size behaviour of the ground state 
and low energy excitations of %e=-K for even L. (h,h) are the predictions from the At (5^) 
minimal model. We have also indicated the D(D%) sector in which the state appears and its 
degeneracy The operator content of the sector irj, is obtained from that of 7r c by interchanging 
h and h. 



D(D 3 ) 


j^num 


(h,h) 


spin 


degeneracy 


TT a © VT b 


0.000000(1) 


(0,0) 





1 x 32- 1 




0.801(3) 







1 x 32- 1 




1.80(1) 


(1 1\ (1 

v5' 5/> V5' 5^ 


±1 


1 x 32- 1 




0.4668(2) 


f_L 2x 


1 

3 


2 x 3^ -1 




0.666666(1) 


(l.o) 


2 
3 


2 x 32- 1 


7T e © 7Tf 


0.13334(6) 


-) 

Vl5' 15 ' 





4 x 3*- 1 




1.33333(3) 


(2 2x 
V3' 3^ 





4 x 32 - 1 



On the other the irreps 7r a and 7Tf, can be considered separately but as we have massive degen- 
eracy at 6 = 7T we have the same exponents appearing in both the 7r and 7T[, sectors, so it is 
convenient to combine them. The degeneracy of each excitation is also listed in the table. These 
degeneracies are determined by the symmetry sector, the number of states in that symmetry 
sector and the pairing multiplicities of those states (which we discuss later). We also classify 
the different symmetry sectors in terms of Bethe root configurations, in particular the number 
of roots at ±00, summarised in the Table [2j Due to the symmetry of the Bethe equations we 
can map every Bethe root configuration to another, {xj} — > {— Xj}. This mapping does not 
change the symmetry sector for states within ir a © or 7r c © Ttf, but maps states from 7r c to 7r 
and vice-verse. Hence the (h, h) (h, h) correspondence between the ir c and 7r a sectors. 



Table 2: The classification of symmetry sectors in terms of finite and infinite Bethe roots for 
chains of even length. 



Symmetry Sector 


N_ + N + + 2N S 


Ti—oo 


ra+oo 


vr a © 7Tb 


C 










C-l 





1 




C-l 


1 





7T C © 7Tf 


C-2 


1 


1 



The excitations appear here for C even coincide with those from the self-dual ferromagnetic 
3-state Potts quantum chain subjected to either periodic of twisted boundary conditions |55j . 
This stems from both models being constructed from the set of solutions to the star-triangle 
equation [16], albeit with different limits applied. The excitations in the 7r © tt^ and ir c © 
7Tj sectors belong to the periodic 3-state Potts chain and are distinguished by translational 
properties, where as the excitations in the ir c and 7r a sectors correspond to the 3-state Potts 
chain with twisted boundary conditions allowing for the non-half-integer spins observed. The 
excitations appearing in the periodic Potts chain have also been determined using Bethe ansatz 
methods by Albertini et al. [2]. 

We have also determined the excitations for chains of odd length. In this case the states can 
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be classified by the eigenvalue of the state under the D3 rotation operator, a. The excitations 
are given Table |3| 



Table 3: As Table [T] for odd C. Symmetry is classified by the action of the D3 rotation a. 



a X^ um (h, h) spin degeneracy 



1 


0.125000(5) 
0.42502(2) 
0.92490(6) 
1.625000(1) 


(o,|) 

(2 M 
v 5 ' 40 1 
( 2 21 x 
v 5 ' 40 1 

(o,f) 
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8 

3 
8 

1 

8 

13 
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C-l 

1 x 3— 

C-l 

1 x 3— 

C-l 

1 x 3^ 

C-l 

1 x 3^r 




0.091665(2) 
0.59168(7) 
0.791667(1) 


-) 

U5 ' 40/ 
(i 21^ 
U5 ' 40 / 
(2 1\ 
V3' 8/ 


1 

24 
11 
24 

13 

24 


£-1 
2 x 3— 

C-l 

2 x 3— 

C-l 

2 x 3~ 



These excitations are not present in the 3-state Potts model since the equivalence to this model 
is restricted to L even. We were able to again classify the symmetry sectors in terms of Bethe 
roots as presented in Table |4j We found that Ti—oo — in the case of odd length chains. 



Table 4: The symmetry sectors classified in terms of finite and infinite Bethe roots for chains 
of odd length. 



Symmetry Sector 


N_ + N + + 2N S 


Tl— 00 


n+oo 


1 


C 








UJ, LU^ 1 


C-l 





1 



4.1.2 Spectrum of Hg=o 

Using analytical and numerical techniques we found that the ground state energy is, 



En 



1 

2^ 



2^/3 3 

— h - 

9 4 



(24) 



Using the Fermi- velocity we find that the central charge is 1, which does not uniquely define a 
conformal field theory. The field content of the theory is obtained from the finite size spectrum. 
One method of determining the the finite size spectrum is using the dressed charge formalism 
(see appendix) leading to the identification of the Z4 parafermion theory |58| [26] coinciding 
with the anti-ferromagnetic 3-state Potts model. The allowed conformal weights for this theory 
are 



h,h G 



1(1 + 2) 
24 



m 

T6 



< m < I < 4, I = m(mod2; 



1 119 3 

' 16' 12' 3' 16' 4' 

Alternately we can solve the Bethe equations directly and determine the scaling behaviour of the 
low-lying excitations. Here we must consider C = 0,1,2,3 (mod 4) separately, see Tables [5]{7 
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below. In particular, we find that the finite size gap of the lowest states for £ = C (mod 4) ^ is 
determined by an (anti-)chiral Zk=A spin field with conformal weight hi = £(k — £)/(2k(k + 2)). 



Table 5: Scaling dimensions Xq extrapolated from the finite size behaviour of the ground state 
and low energy excitations of Hq for C = (mod 4) (the error of the extrapolation is smaller 
than the last displayed digit). (h,h) are the predictions from the Z4 parafermionic CFT. For 



D(D 3 ) 


ynum 


(h,h) 


spin 


degeneracy 




0.000000 


(0,0) 





C 1 
1 X 32" 1 




0.333332 


(o,|) 


1 

3 


2 x 


VT e © 7Tf 


0.166667 
0.666667 


-) 

U2» 12/ 
V3' 3/ 






4 x 3^ _1 
4 x 3#-! 



Table 6: As Table [5] for C = 1 (mod 4). The excitations for chains with length C = 3 (mod 4) 
have the same exponents but the opposite spin. Symmetry is classified by the action of the -D3 
rotation a. 
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spin 


degeneracy 
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0.062500 
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16 


C-l 

1 x 3— 
c-i 
1 x 3— 




0.562500 


(&°) 


9 

16 
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(J, 

U6' 4/ 


11 
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C-l 

1 x 3^ 




0.145833 


-) 
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C-l 

2 x 3— 




0.395833 


U6' 3/ 


13 
48 


C-l 

2 x 3— 

C-l 

2 x 3~ 




0.645833 


(± X) 
U6> 12/ 


23 
48 



Table 7: As Table [B] but for C = 2 (mod 4). 



D(D 3 ) 


ynum 
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spin 


degeneracy 


TT a © 7Tb 


0.750000 


(0,|) x 
2, (f,0) x 2 


=c 4 


C 1 

1 x 32" 1 




0.083333 
1.083333 


(0,^) 

(3 In 
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1 

12 
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2 x 3^ _1 
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7T e © 7Tf 
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V 12 ' 3/' V3' 12J 




4 x 3*- 1 



As with the previous case we again can partition the excitations in into residual symmetry 
sectors. These sectors are still characterised by the number of Bethe roots at ±00 as described 
in Table H 

Comparing these excitation to the anti- ferromagnetic 3-state Potts chain [2,32] (for C even), 
we find that the excitations in the n a © 7T& and 7r c © irt were previously identified and restricted 
to the n+oo = n_oo case. We were unable to find any literature dealing with this end of the 
twisted 3-Potts model, however, we expect that the excitations in that case to match those 
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appearing in the 7r c and 7r 5 sectors. Similarly the excitations for odd C have not previously 
been studied. 



4.1.3 Pairing rules and discussion 

Through studying the two above cases of 6 = 0, tt we are able to deduce the behaviour for all 
the level crossings, i.e. when 9 is a multiple of |, as the energy spectra of and H.^ are 
identical, although the momenta are of opposite sign. The analysis of the model is extended 
to generic 6 by determining pairing rules. Since the symmetry sectors identified above refer to 
the complete system two levels of can pair if and only i^jthey belong to the same 

symmetry sectors^] |21| . Pairing is uniform within symmetry sectors (shown in Table ^ and 
can be used along with Equations ( 9|11|23(24 ) to determine the complete energy. The pairing 



of levels also ensures the locality of the physical fields by restricting spins to either whole-, half- 
or quarter- integer. For example, for £ = (mod 4) with = we find the excitations (^, |) 
from Tin and (0, ^) from H% (same spectrum as %$) both belong to the 7r c symmetry sector, 
thus they pair together twice and have a spin of ^ — = 



Table 8: The pairing multiplicities, m p , of the periodic spin chain for any two energies of 7£W 
and that belong to the same symmetry sector. If the energies do not belong to the same 
symmetry sector then they do not pair, i.e. the pairing multiplicity is zero. 







c = 


(mod 2) 




C = 1 (mod 2) 


Sector 








7T E © 7Tf 


1 w^w^ 1 


m p 


1 


2 


2 


4 


1 2 



The complete model can be summarised by Figure [3} 
4.2 Periodic Fusion Path Chain 

By construction this chain must share bulk properties with the periodic spin chain of the 
previous section, including energy per unit lattice site, Fermi velocity and central charge. Fur- 
thermore we find that many of the energies appearing in the periodic fusion path chain are 
also present in the periodic spin chain. Specifically, these are the energies in the 7r a © 7T(, and 
7r e © 7Tf symmetry sectors and are numerically exaclQ However the degeneracies of the identical 
energies in the two formalisms do not match. Amongst the energies appearing in both are the 
the ground states, when realised, given by Equations ( p3"| ) and (24) for Ti^ and %q, respectively. 



4.2.1 Spectrum of l~Le=-K 

We are still in the M( 56 ) universality class (c = |). Using numerical methods we find the 
excitations as given in Table [9} 

More accurately, the only if direction is true because of the symmetry sectors being a property of the 
complete model while the if direction has been observed to be true. 

6 If an excitation (h, h) of iiH^ 1 belongs to a certain symmetry sector then the corresponding excitation (h, h) 
of will belong to the same symmetry sector. We will later see this is not the case for the braided and open 

chains. 

7 At least up to system sizes of C — 10. Systems with longer chain lengths were not numerically computed. 
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Figure 3: Conformal field theory description of the periodic integrable D{D%) symmetric model 
for generic 9. Level crossing are described by minimal models, either Z4 parafermions or Air^^, 
while the regions in between are described by a product of the theories of the two adjacent level 
crossings. Therefore regions I and IV correspond to a Z4 <g> -M( 5 6 ) theory, while regions II and 
III correspond respectively to eg) Z4 and -M( 5j6 ) <S> -M(5,e) theories. 



cos(#) 




* sin(0) 



For energy levels appearing in both formalisms we have used the numerical Bethe ansatz results 
from the previous section to determine the exponent, while for the new energy levels we have 
used data from the numerical diagonalisation of the Hamiltonian for chains sizes up to 10 sites. 
The new excitations appearing in this chain also do not correspond to excitations of the 3-state 
Potts chains subject to the boundary conditions previously studied [S3]. Moreover because the 
formulation of this chain is reliant on D{D^) symmetry, which is not present in the usual 3-state 
Potts local Hamiltonian, it is reasonable to expect that these excitations won't appear for any 
formulation of the 3-state Potts model. The periodicity of the chain also breaks the D(D^) 
invariance of the model |19| and we are unable assign symmetry sectors to the excitations as 
we did before. 



4.2.2 Spectrum of H e =o 

Using numerical methods we can find the excitations of the occurring for the Hamiltonian %e=o- 
This is done for even L and split according to if ^ is even or odd. The excitations are, 

( h h) G J {(0.0),(A,A)*,(A,^),(i, j)}, £^0mod4, 
1 ' J I {(is, (0, !), (|,0), (i, I)}, C = 2mod4, 

where we have again have used a * to denote new excitations. As the spectra need to be 
considered differently for i odd and even, we have effectively split our number of data points 
for each excitation. Thus the accuracy here is lower than the excitations for Tie^. 



4.2.3 Pairing rules and discussion 

In contrast to the periodic spin formulation the pairing rules for the periodic fusion path chain 
are not of a simple form. Pairing is no longer transitive implying the pairing table would have 
no block structure. Additionally the symmetry sectors used to describe the periodic spin chain 
can not be used to describe this case as the periodic closure of the fusion path chain breaks 
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Table 9: Scaling dimensions X v extrapolated from the finite size behaviour of the ground state 
and low energy excitations of H n for C = (mod 2). The values with a * correspond to 
exponents not appearing the periodic spin formulation and therefore also not in the 3-state 
Potts model. Scaling dimensions marked with a i" have no estimation of uncertainty due to 
a lack of data points. Each combination of (h, h) appearing in the second column has the 
degeneracy given in the last column. 
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the D(D%) symmetry in the spin sense. However, for the anyon model it is possible to define 
a conserved topological charged bases of the D(D^) F-moves aforementioned [TTJ, [27]. The 
different topological sectors are labelled by the irreps in the decomposition of tt® £ . A better 
understanding of these topological symmetry sectors may lead to insight into the pairing rules. 

The identification of the excitations has utilised techniques which numerically diagonalise 
the Hamiltonian. This is clearly limited by computational power. As the periodic fusion path 
chain shares energies with the periodic spin chain we have confidence in the ones we have 
successfully identified. However, the excitations which are do not appear in the periodic spin 
chain are difficult to determine within any great accuracy. To ensure that we have identified 
all excitations of the level crossings and answer the open problems associated with this model, 
such as pairing rules, more advanced techniques would be required such as using integr ability 
reliant techniques associated with RSOS models or utilising better numerical methods for the 
diagonalisation of Hamiltonian. 



4.3 Braided Chain 
4.3.1 Spectrum of T-Le=-K 

The full spectrum of the Tle=-K braided chain is a subset of the spectrum of the He=n periodic 
spin chain. In particular we find that the energies present are those that appeared in the 
symmetry sectors ir a © 7r& and ir c of the periodic chain for even C and in the a = 1 for odd 
C This corresponds to Bethe root configurations with n_oo = 0. Therefore the low- lying 
excitations for this chain correspond to the pairs of conformal weights, 

(hh) e / {(0,0), (§, §), (§, |), (|, |), |), (§,0)} , C even, 
( ' } ' I {(0,|), (1,^,(1, |), (0,|)}, C odd. 

The braided chain has the full D(D^) symmetry implying we could classify excitations in terms 
of the irreps, however, we find that because of the difference in global symmetry we now have 
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energies appearing across multiple symmetry sectors. We also remark that since there are 
less excitations each of them has a higher degeneracy than the same excitation in the periodic 
chain. For chains of even length the degeneracies are now either 3x3^ _1 or 6 x 3^ _1 depending 
on which symmetry sectors they are in. While for chains of odd length every excitation has 

C-l 

degeneracy 3 x 3 2 . 
4.3.2 Spectrum of Uq=q 

As %e = o = —T-Le=ir w e again have the excitations appearing also were identified in the periodic 
section. The ones that appear are, 



(h,h) e < 



{(0,0), (0±)}, £ = 0mod4, 

{(0,f),(|,0),(0,i),(|,i)}, £ = lmod4, 

{(1 0),(&,0),(i I)}, £^2mod4, 

I {(O4), (l,0),(i,0), (!,§)}, £ = 3mod4. 

4.3.3 Pairing rules and discussion 

Like the periodic case by studying the braided chain at 9 = 0, ir we can extrapolate finite size 
data to determine the behaviour at all of the level crossings. We observe important differences 
between the periodic and braided models. We find that an energy of Ti^ will appear in multiple 
symmetry sectors and similarly for the energies of%( 2 ). However, we observe that for a specific 
exciation, (h, h), the symmetry sectors it appears in for 9 = 0, tt are not necessarily the same as 
the ones it appears in for 9 = §, 4f- Excitations appear in multiple symmetry sectors in such 
a way that every two levels pair. For ir < 9 < 4? and £ even we construct the pairing table 



(Table 10) of the lowest excitations. In this table we can clearly see excitations of Hq=k and 

Ha—in. appearing in multiple symmetry sectors in such a way that every excitation pairs with 
°~ 2 

all others. We should also remark that we could have formulated this in the fusion path basis 
which would have provided an equivalent model up to degeneracies. 



Table 10: Pairing table for the lowest excitations of the periodic spin chain for tt < 9 < 4?. 
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4.4 Open Chain 

We now consider the open chain. At level crossings conformal invariance predicts that the 
lowest energies of the open chain with free ends, i.e. xf = X2 = 0, have the scaling behaviour 
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given by 



E = e 00 C + h + 7n ^(-^ + h + n) +o(C~ 1 ), 



where c is the central charge, vf is the Fermi- velocity, h is a conformal weight and n is a 
non-negative integer. 

4.4.1 Spectrum of %e= v 

The ground state (h = 0) energy is, 
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As is the case with the periodic chain we have that vp = 3 and c = |, yielding the same CFT 
as expected. From our numerical solutions to the Bethe ansatz equations we found the lowest 
energies corresponded to conformal weights 



h G 



£ even, 



{|,f}, £odd. 



and descendants thereof. 



4.4.2 Spectrum of ^6»=o 

The ground state of the open chain for 6 = is, 
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which gives up = | and c = 1, in agreement with the periodic chain. Numerical data showed 
that the lowest excitations appearing correspond to the conformal weights 



(h,h) G 



along with some their descendants. 
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I 16' 16/ ' 
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£ = 0mod4, 
£ = lmod4, 
£ = 2mod4, 
£ = 3mod4. 



4.4.3 Pairing rules and discussion 

This model, like the braided version, has the full global symmetry of the algebra D{D%) and 
also has all energies pairing. 



5 Discussion 

In this paper we have analysed the low energy spectrum of the integrable D{D^) symmetric 
chain subject to various boundary conditions. In the spin chain formulation the Hamiltonian 
derives from a commuting two-parameter transfer matrix of a vertex model and the eigenvalues 
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can be obtained by Bethe ansatz methods. We have constructed a related class of models with 
local D(Ds) symmetry using the fusion path formulation: these models, too, are integrable as 
they can be obtained from a solution to the Yang-Baxter equation for a face (or RSOS) model. 
For open and braided boundary conditions the two formulations of the model are equivalent. 
For periodic closure, however, the fusion path chain differs from the spin chain by boundary 
terms. 

Based on our analysis of the spectrum of these models we have identified the conformal 
field theories providing an effective description of the low energy modes to contain two sectors 
from - depending on the parameter 6 - the minimal model At (5^) (the three-state Potts model) 
and the Z4 parafermion. The physical fields are products of operators from these sectors. The 
individual factors can carry fractional (non- integer or non-(para)fermionic) spins implying the 
appearance of Virasoro characters in the partition function of the model which have not been 
discussed in the context of -M( 5 6 ) or Z4 alone. 

The locality of physical fields in the model is guaranteed by pairing rules. This situation is 
similar to other models with several gapless modes propagating with different Fermi velocities 
[H [2"5t 124] . We have to emphasise, however, that in the present models this factorisation 
of different modes is exact already for finite chains and on all energy scales, unlike in say the 
separation of spin and charge degrees of freedom observed within the low energy spectrum of the 
one-dimensional Hubbard model where the coupling between the sectors becomes manifest in 
subleading corrections to scaling. Another difference for the model studied here is that the two 
sectors of the effective theory are not related to subalgebras of the global symmetry of the model. 
Therefore, to establish the pairing rules and corresponding multiplicities we have resorted to 
numerical studies of small systems together with counting arguments for the total number of 
states of the system. For the spin chain formulation we found that the pairing is determined 
the boundary conditions and can be related to the residual symmetry of a given eigenstate. 
It is also reflected by the appearance of infinite roots appearing in the configurations solving 
the Bethe equations of the model. The spectrum of the periodic D(D$) model in the fusion 
path formulation also shows pairing on all energy scales. Unlike in the spin chain formulation, 
however, the pairing is not transitive and cannot be described based on residual symmetries 
as in the spin chain. The determining of the pairing rules and as to whether these rules are 
connected to topological invariants involving the D(D^) F- moves is an interesting open problem 
of this model. 

As a first step, numerical methods could be used to form conjectures. Ultimately, however, 
it would be desirable to obtain the complete picture starting from the integrable structures un- 
derlying this model. Apart from the understanding of pairing rules this concerns in particular 
the solution of the spectral problem of the anyon chain: the fact that the two formulations of 
the periodic chain have exponentially (in system size) many energies in common suggests that 
there should exist a formulation of the eigenvalue problem in terms of Bethe type equations for 
the fusion path chain. We shall address this problem in future work. 
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A Dressed charge formalism 

Following [61 [52j [2H US] the finite size energy gaps of rle=o for periodic boundary conditions 
are 

AE(AN±, Q±) = 2 ^v F (](AN + ,AN^)(E r E)- 1 (AN + ,AN^ 

C U v (25) 

+ (Q+, Q_)(~ T ~)(Q+, Q_)+M)+ o(C~ l ) . 



(Af being a non- negative integer). Taking into account that only i-strings are present in the 
ground state the 2x2 dressed charge matrix H = £(x)\ x=00 is obtained from the linear integral 
equation 



k(x, |) k(x, s ) J ' ' cosh(x) - cos (2-irt)' 

Using Wiener Hopf techniques the dressed charge matrix can be expressed in terms of the 
Fourier transform of the kernel matrix giving 



(l -£(« = (>)) =( [ f . (27) 
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Hence the scaling dimensions and conformal spins of primary operators in the effective field the- 
ory for 710=0 m terms of the quantum numbers AN± and Q± characterising the corresponding 



excitation (22) are 

X = i {(AN + f - AiV+AAL + (AiV_) 2 ) + ((Q + ) 2 + Q + Q_ + (Q_) 2 ) , 
s = ~(Q+AN + + Q_AN_) . 



(28) 



The AN± = N± — ^ ± j correspond to the change in number of i-strings (subject to the 
condition that the total number of roots is C) as compared to the thermodynamic ground state 
while the Q± determine the momentum of the excitation. For a configuration with n± Bethe 
roots at ±oo they can take discrete values 

Q± = ^{n+oo ~ "--oo) (mod 1). 

For a given solution {x^} of the Bethe equations the Q± can also be determined numerically 
using the counting functions defined above: 

1 N± 
± k=l 

Due to the discrete set of possible values for Q± the data for systems with C < 10 are sufficient 



to identify the quantum numbers for the lowest finite size gaps of 1-Le=o, see in Table 11 The 



observed dimensions support our identification of the critical theory with a Z4 parafermion. 



Table 11: The lowest excitations of T-Lq in terms of the quantities AN± and Q± for different 



chain lengths. 
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